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20  ABSTRACT  (Continue  on  reverse  side  it  necessary  end  identify  by  block  number) 


1  Theoretical  and  experimental  studies  of  planar  detonation  waves  in 
Amatex  20  were  performed  to  obtain  a  better  understanding  of  the  role  of 
NHijNGj^-^AN)  in  the  detonation  process  in  explosives  formulated  from  RDX,  TNT, 
and  AN.  1  >  a  u 

The  theoretical  work  was  performed  to  provide  a  more  definitive  and 
realistic  treatment  of  steady-state  detonation  waves  in  such  composite 
explosives.  The  treatment  developed  for  such  waves  extends  the  classical 
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20.  ABSTRACT  (continued) 

-  ^ Zelodovich-von  Neumann-Doering  (ZND)  model  for  ideal  detonation  into  a 
prototype  model  for  nonideal  detonation.  The  kinetic  processes  respon¬ 
sible  for  nonideal  behavior  were  assumed  to  be  incomplete  decomposition 
of  the  slowest  reacting  explosive  component,  and  the  incomplete  attain¬ 
ment  of  chemical  equilibrium  among  the  detonation  products  from  the 
different  explosive  components.)  The  constitutive  relationship  for  the 
shocked  reacting  explosive  was constructed  with  different  equations  of 
state  for  the  explosives  and  their  products  and  enough  reaction  coor¬ 
dinates  to  account  for  this  type  of  nonideality.  This  constitutive 
relationship  was  combined  with  the  Rankine-Hugoniot  (RH)  relationships, 
and  the  Chapman-Jouguet  (CJ)  sonic  condition  to  derive  equations  for 
(1)  the  CJ  parameters  in  and  the  propagation  velocity  of  a  steady-state 
nonideal  detonation  wave,  and  (2)  the  particle  velocity  in  terms  of  the 
reaction  coordinates  in  a  steady-state  nonideal  detonation  wave.  Detona¬ 
tion  parameters  and  Lagrange  particle  velocity  histories  were  calculated 
for  nonideal  detonation  waves  in  an  explosive  modeling  Amatex  20. 

/The  experimental  work  was  performed  to  obtain  a  more  complete  charac¬ 
terization  of  the  detonation  process  in  Amatex  20.  The  mean  particle  size 
of  AN  in  the  Amatex  20  charges  was  increased  from  the  baseline  value  of 
-0.7  mm  to  a  value  of  -2  mm.  Multiple  Lagrange  particle  velocity  gage 
experiments  were  performed  and  analyzed  to  obtain  a  quantitative  descrip¬ 
tion  of  planar  detonation  in  this  coarse-grained  Amatex  20. *•  The  experi¬ 
mentally  determined  CJ  parameters  for  this  coarse-grained  Amatex  20  _ 

(D  =  6.85  mm/ys,  u^  =  1.67  mm/ys,  pCJ  =  182  kbar,  v^  =  0.4742  cm3/g, 

KCJ  =  3.10,  TCJ  ~  0.63  ys)  and  those  determined  for  the  baseline  Amatex  20 

(D  =  6.95  mm/ys,  uCj  =  1.94  mm/ys,  p^  =  217  kbar,  vCj  =  0.4470  cm3/g, 

Kr  =  2.58,  x  -  0.5  ys)  clearly  demonstrate  that  the  particle  size  of  AN 

influences  the  detonation  process  in  Amatex  20.  Kamlet's  semi-empirical 
method  for  calculating  detonation  parameters  was  modified  and  used  to 
interpret  the  differences  between  these  sets  of  CJ  parameters.  Application 
of  this  method  led  us  to  formulate  the  following  hypotheses  about  the 
nonideal  planar  detonation  process  in  Amatex  20:  (1)  AN  with  a  mean 

particle  size  of  0.7  mm  decomposes  completely  in  the  reaction  zone,  but 
the  decomposition  products  of  RDX,  TNT,  and  AN  do  not  attain  chemical 
equilibrium,  and  (2)  AN  with  a  mean  particle  size  of  -2  mm  does  not 
decompose  completely  in  the  reaction  zone,  and  the  detonation  products 
of  RDX,  TNT,  and  AN  do  not  attain  chemical  equilibrium. 
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I  INTRODUCTION 


The  overall  objective  of  the  research  was  to  obtain  a  basic 
understanding  of  ammonium  nitrate  (AN)  In  the  planar  detonation  process 
In  composite  military  explosives  formulated  from  RDX,  TNT,  and  AN. 
Theoretical  and  experimental  studies  were  performed  to  extend  the 
previous  work^  performed  under  Contract  No.  DAAG29-76-C-0033. 

The  theoretical  work  was  performed  to  provide  a  more  definitive  and 
realistic  treatment  of  detonation  in  composite  explosives  such  as  Amatex 
20.  We  addressed  the  general  problem  of  detonation  in  a  composite 

explosive  X  fabricated  from  three  explosive  components,  X].,  X2,  and  X3* 
It  is  necessary  in  approaching  this  problem  to  identify  and  model  rate 
processes  that  may  influence  the  detonation  process  and  give  rise  to 
nonideal  behavior.  The  model  for  planar  detonation  in  such  composite 
explosives  constructed  in  this  report  is  based  on  the  premise  that 
certain  chemical  reactions  do  not  proceed  to  completion  in  the  steady- 
state  reaction  zone.  In  this  case,  nonideal  behavior  arises  because  not 
all  the  available  chemical  energy  is  used  to  support  the  propagation  of 
the  detonation  wave. 

The  experimental  work  was  performed  to  obtain  a  more  complete 
characterization  of  the  detonation  process  in  Amatex  20.  The  mean 
particle  size  of  the  AN  in  the  Amatex  20  charges  was  increased  from  the 
baseline  value  of  ubout  0.7  mm  to  a  value  of  about  2  mm.  Multiple 
Lagrange  particle  velocity  gage  experiments  were  performed  and  analyzed 
to  obtain  a  quantitative  description  of  planar  detonation  in  Amatex  20 
containing  this  coarse-grained  AN.  Significant  differences  in  the 
values  of  the  detonation  parameters  in  the  charges  containing  base-line 
and  coarse-grained  AN  clearly  demonstrate  that  the  detonation  process  in 
Amatex  20  is  influenced  by  the  particle  size  of  the  AN.  The  method 
developed  by  Kamlet?-^  to  calculate  CJ  parameters  was  modified^  and  used 
to  interpret  these  differences. 


II  THEORETICAL  STUDY 


Our  theoretical  study  is  based  on  the  assumption  that  planar  one¬ 
dimensional  detonation  in  our  composite  explosive  can  be  treated  as  a 
Zeldovich-von  Neumann-Doering^-®  (ZND)  wave.  In  this  case,  states  in  a 
steady-state  detonation  are  governed  by  the  Rankine-Hugoniot  (RH)  jump 
conditions^  expressing  the  balance  of  mass,  momentum,  and  energy  in  a 
steady-state  flow.  The  flow  becomes  sonic  at  the  Chapman-Jouguet  (CJ) 
point  where  the  global  energy  release  rate  (GERR)  becomes  zero,  and  when 
all  the  reactions  contributing  to  the  (GERR)  are  exothermic,  the  CJ 
point  is  at  the  end  of  the  reaction  zone.  Exothermic  reactions 
considered  to  contribute  to  the  GERR  are  as  follows: 

•  The  global  decompositions  of  the  explosive  components 

^1 ,  X2,  and  X3  into  their  detonation  products  in  local 
chemical  equilibrium. 

•  Reactions  among  the  different  detonation  products  of  X^ t 

x2»  and  X3  as  they  mix  and  react  to  approach  global 
chemical  equilibrium. 

We  are  now  in  a  position  to  define  "ideal”  and  "nonideal”  detonation. 

A  steady-state  ZND  wave  is  said  to  be  ideal  when  all  the  global 
decompose  ..on  reactions  of  the  explosive  components  proceed  to 
completion  and  the  products  of  these  reactions  attain  global  chemical 
equilibrium.  Otherwise  it  is  said  to  be  nonideal.  This  definition  is 
compatible  with  the  one  we  used  previously^-  and  leads  us  to  assume  that 
the  two  major  factors  responsible  for  nonideal  behavior  of  our  composite 
explosive  are  as  follows: 

(1)  The  incomplete  decomposition  of  a  slower  reacting 
explosive  component  such  as  AN. 

(2)  The  unattainment  of  chemical  equilibrium  among  the 
detonation  products  from  different  explosive  components. 
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We  envisage  the  first  factor  to  be  important  when  the  decomposition  of 
the  slowest  reacting  component  is  not  self-propagating  and  terminates 
when  the  decomposition  of  the  faster  reacting  components  has  gone  to 
completion.  We  envisage  the  second  factor  to  be  important  when  the  time 
scale  for  the  detonation  products  to  mix  and  react  exceeds  the  time 
scales  for  the  decomposition  reactions.  The  constitutive  relationship 
for  reacting  composite  explosive  used  in  a  treatment  of  nonideal  detona¬ 
tion  must  include  enough  reaction  coordinates  to  account  for  nonideality 
arising  from  factors  (1)  and  (2). 

The  theoretical  model  for  nonideal  detonation  presented  in  this 
report  was  formulated  to  provide  a  better  understanding  and  a  means  for 
approaching  the  problem  of  steady-state  detonation  in  composite  explo¬ 
sives.  With  Xlf  X2,  and  X3  considered,  respectively,  to  be  like  RDX, 
TNT,  and  AN,  it  was  necessary  to  construct  a  constitutive  relationship 
for  the  reacting  explosive  with  the  four  reaction  coordinates  X^ , 

X^,  and  X^.  The  first  three  coordinates  account  respectively  for  the 
global  decompositions  of  X^ >  x2,  and  X3,  and  the  fourth  accounts  for  the 
recombination  of  the  oxygen  produced  by  X3  and  the  carbon  produced  by  X3 
and  X2.  This  constitutive  relationship,  the  RH  jump  conditions,  and  the 
CJ  sonic  condition  were  combined  to  generate  the  equations  governing  the 
steady  reaction  zone  and  thereby  demonstrate  the  dependence  of  the  deto¬ 
nation  process  in  composite  explosives  on  incomplete  reaction  and  on  the 
equation  of  state  of  the  condensed  explosive.  The  expression  relating 
the  particle  velocity  and  the  reaction  coordinates  in  a  steady-state 
reaction  zone  was  also  derived.  This  expression  ias  combined  with  model 
reaction  rate  expressions  to  generate  Lagrange  particle  velocity 
histories  and  to  investigate  the  dependence  of  particle  velocity  on  the 
relative  rates  of  the  exothermic  reactions. 

Governing  Equations 

Flow  Equations 


We  first  introduce  the  notation  used  to  describe  our  one-dimensional 
flow.  Time  is  denoted  by  t,  Lagrange  distance  by  h,  specific  volume  by  v 


particle  velocity  by  u,  pressure  by  p,  sound  speed  by  c,  and  specific 
energy  by  e;  the  subscript  x  denotes  explosive,  and  the  superscript  o 
denotes  the  constant  state  ahead  of  the  wave.  Neglecting  the  initial 
pressure,  the  flow  in  our  steady-state  ZND  wave  is  governed  by  the  RH 
jump  conditions  written  as 


Vx  0>  - 

u) 

(1) 

Du/v°  = 

(D  -  u)u/v 

(2) 

1 

f  o  ■>  12 

(3) 

=  2  P 

x  '  2 

where  D  denotes  the  propagation  velocity.  We  use  H  as  a  subscript  or  as 
a  superscript  to  denote  the  shocked  state  at  the  wave  front,  and  we  use 
CJ  as  a  subscript  to  denote  the  CJ  point.  Then  the  equations  for  the 
shocked  state  are  obtained  by  setting  v  ■  vj|,  u  =  uH,  p  =  pjj,  and  e  =  e 
in  Eqs.  (1)~(3),  and  the  equations  for  the  CJ  state  are  obtained  by 
writing  the  CJ  condition  as 


D  "  UCJ  “  cCJ 


(4) 


and  setting  v  =  v^j,  u  =  u^j,  p  =  p^j,  and  e  ■  e^j  in  Eqs.  (l)-(3).  The 
adiabatic  inviscid  flow  in  the  wave  behind  the  shock  is  governed  by  the 
following  expressions  for  the  laws  of  conservation  of  mass,  momentum, 
and  energy: 


(5) 

(6) 
(7) 


The  flow  and  chemistry  in  our  ZND  wave  are  coupled  by  the  constitutive 
relationship  relating  e  and  the  reaction  coordinates  X^ ,  X^,  X ^  and  \^, 
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xas 


Equations  for  the  Explosive 

Following  the  papers  by  Kamlet  et  al.  ,  q  we  write  the  molecular 

formulas  of  our  explosive  components  X^,  X2,  and  X3,  respectively,  as 

C  H,  N  0.  with  i  =  1,  2,  and  H,  N  0,  .  We  also  assume  that  the 
ai  bi  ci  di  b3  c3  d3 

decompositions  of  X^  and  X2  satisfy  the  H2O  -  CO2  arbitrary  equilibrium 

condition.  The  H2O  -  CO2  arbitrary  represents  N2,  H2O,  and  CO2  as  being 
the  only  important  detonation  products,  with  1^0  having  priority  in 
formation  over  CC^.  We  accordingly  write  the  decomposition  of  X^,  with 
i  =  1,  2,  as 


C 

a 


i 


(8) 


Extending  the  notion,  we  assume  that  the  decomposition  of  X3  satisfies 
the  H2O  -  O2  arbitrary  so  that  N2,  H2O,  and  C>2  are  the  only  detonation 
products,  with  H2O  having  priority  in  formation  over  C>2«  We  accordingly 
write  the  decomposition  of  X3  as 


c3 

~2~  N2 


H  20  + 


(9) 


The  reaction  coordinates  X.^,  and  ^3  are  used  to  describe  the  decom¬ 
positions  specified  by  Eqs.  (8)  and  (9),  and  the  reaction  coordinate  \ ^ 
is  used  to  describe  the  reaction 


d3  b3  \ 


CO , 


(10) 


with 


(d 


3/2  "  b3/4 


)  <  ,S7 

1  -  i  =1 


<*1  ‘  di/2 


+  bj/4) 


that  may  occur  among  the  decomposition  products  of  X^ f  X2,  and  X3. 
Because  of  the  H20  -  CO2  and  H2O  -  02  arbitrary  assumptions,  only  one 
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reaction  coordinate  is  needed  to  account  for  subsequent  reactions  among 
detonation  products.  With  the  standard  heat  of  formation  of  a  species 
Sj  denoted  by  (AH°)g  ,  the  standard  heats,  Q^,  liberated  in  our 
exothermic  reactions  Eqs.  (8)-(10)  are  given  by  the  equations 

«i  ■  -<V2>  (s“?>h20  -  <di/2  '  V4>  <“?>C02  +  (“?)Xl  <U) 

with  i  =  1,  2 

q3  .  -(bj/2)  <™0f)x3 

%  '  -<d3/2  -  b3/4><al°f>C02 

where  (AH°)  *  -57.8  kcal/mole  and  (AH°)_  *  -94  kcal/mole.  We  now 

t  t  CO 2 

introduce  the  notation  used  to  describe  our  explosive  mixture  and  its 

reaction  products.  We  let  the  subscript  p  denote  detonation  products, 

let  Mg  denote  the  molecular  weight  of  a  product  species  Sj,  and  we 

let  v°  ,  e°  ,  N^,  and  denote  the  initial  specific  volume,  the  initial 

specific  internal  energy,  the  number  of  moles,  and  the  mass  fraction  of 
the  explosive  in  our  composite  explosive  X.  Then  and  are 
related  by  the  equation 


(12) 

(13) 


N .  -  ct./M  i  -  1,  ...  3 

i  i  Xj 


(14) 


With  ,rr  N,M  *1.  We  assume  that  the  mixture  of  explosives  is  ideal 
1*1  i  x^ 

and  write  the  following  equations  for  its  specific  volumes,  specific 
internal  energy,  and  specific  heat  of  formation: 


v 


o 

x 


e 


o 

X 


V  3  o 
i*n  ai  V 


(15) 

(16) 
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where  (Ah”)  *  (AH?)  /M0  .  If  we  assume  that  the  initial  pressure  is 

Sj  Sj  Sj 

zero,  Eqs.  (16)  and  (17)  are  identical  because  (Ah?)  ■  e°  .  We  let 

f'xt  x^ 

p^  denote  the  mass  fraction  of  a  decomposition  product,  Sx,  of  X^  when 

j 

has  decomposed  completely.  Then  under  the  H2O  -  CO2  and  H2O 
arbitrary  assumptions  it  follows  from  Eqs.  (8)  and  (9)  that 

(18) 

and 

+  ^20  +  ^  1 

(19) 

where 

CM 

CD. 

-  Ci"n2/2“k1  1  -  2.  3 

^2o 

’  blMH20/2Mx1  1  ‘  2.  2.  3i 

4o2 

•  (di  -  V2>  MC02/2Mx1  1  ■  3'  2! 

-  [2.J  -  <d1  -  b±/2) ]  Mc/2Mx  1  *  1,  2; 

-  <d3  *  b3/2>  \'2\ 

We  can  then  write  the  standard  heats  q^,  q2,  and  q3  liberated  by  the 
decomposition  of  1  gram  of  X^,  X2 ,  and  X3  as 

A 

-  el  n  (Ah°)u  n  -  pL  (Ah?)  +  (Ah? )  i  -  1,  2 

(20) 

(22) 


«3  ■  -  %20  <a>;>H2 0  +  K>*3 

A 

and  the  standard  heat  q^  liberated  by  the  decomposition  of  X^,  X2>  and 
X3  in  1  gram  of  X  as 


qd 


ai  <*i 


(22) 


We  now  consider  the  reaction  among  the  detonation  products,  Eq.  (10), 
and  write  the  standard  heat  of  this  reaction  for  1  gram  of  X3  as 


^4 


(23) 


where  (3^  ■  (d^/2  -  b^/4)  Thus,  when  all  the  exothermic 

reactions  proceed  to  completion  in  our  composite  explosive,  we  can  write 
the  standard  heat  of  these  reactions  for  1  gram  of  X  as 


q  -  *i+  % 

A  A  A  A 

where  ql  =■  a^,  q2  -  a^,  q3  =  a3q3,  and  q4  =  a3  q4- 


(24) 


Constitutive  Relationship  for  Reacting  Composite  Explosive 

Our  constitutive  relationship  for  reacting  composite  explosive  is 
an  extension  of  the  constitutive  relationship  formulated  previously*-® 
for  an  explosive  with  one  global  decomposition  reaction.  We  treat  the 
explosive  components  and  their  reaction  products  as  a  mixture  of  phases 
that  attains  mechanical  but  not  thermal  equilibrium.  It  is  also  assumed 
that  no  appreciable  amount  of  heat  is  transferred  into  the  composite 
explosive  as  the  decomposition  reactions  propagate  and  consume  the 
explosive  components.  In  this  case,  the  explosive  components  are  com¬ 
pressed  or  released  isentropically  as  the  reactions  proceed,  and  the 
pressure  increases  or  decreases  along  a  particle  path. 


We  assume  for  convenience  that  the  explosives  and  their  products 
are  polytropic  materials,  and  we  assume  for  tractability  that  the  poly¬ 
tropic  indices  of  the  explosive  components  are  equal,  i.e.,  ^  = 

=  Kx,  and  that  the  polytropic  index  of  carbon  is  the  same  as  that  of 
the  other  reaction  products,  i.e.,  =  Kp.  In  contrast  to  classical 

treatments  of  detonation,  however,  we  do  not  assume  that  K  =  K  =  K  so 

x  p 

that  we  can  investigate  the  influence  of  the  equation  of  state  of  the 
condensed  explosive  on  the  detonation  process.  We  accordingly  write  the 
(e-p-v)  relationships  for  the  explosive  components  and  their  reaction 
products  as 


(Ahf0)Xl  + 


(K  -  1) 
Xi 


i  =  1,  2,  3 


=  (Ah”)„  + 


(Kp  -  1) 


We  denote  the  specific  entropy  by  s  and  write  the  mixture  rules  for  the 
reacting  explosive  mixture  as 


e  =  Se  +  e 
x  p 


v  =  Ev  +  v 
x  p 


s  =*  Ts  +  s 
x  p 


z\  '  A  “i(1  ‘  V  v 


sx  *  A  “i(1 '  V  v 


<<>x  -  ViK>x 


K 


(30) 


and 


ep  "  i-1  Xiai0H2O(Ah  f^O  +  i^l  XiaiPC02(/!hf)C02 

+  +  X3X4“3PCO,(Ahf  ^O,  +  (K"  -_T)  <31> 

^  4  Z  Z  p 

The  combination  of  Eqs.  (27),  (30),  and  (31)  after  some  manipulation 
gives  the  constitutive  relationship  for  our  reacting  mixture  as 

e  -  (Ahf  )x  -  ^  \lq1  -  *4*394  +  K  P-  1  "  K?  -  1  (32) 

P  P 

A  A 

with  q±  *  a1q1,  q4  -  a^,  and  =  (Kx  -  Kp)/(Kx  -  1).  The  (e-p-v) 
relationship  for  our  composite  poly tropic  explosive  is  obtained  from  Eq. 
(32)  by  setting  X^  *  X^  “  X^  =*  0,  v  *  vx>  and  e  *  e  ,  and  the  (e-p-v) 
relationship  for  our  polytropic  decomposition  products  is  obtained  from 
Eq.  (32)  by  setting  X^  *  X^  ■  X^  =  1,  v  *  vp ,  and  e  *  ep*  The  classical 
(e-p-v-X^)  equation  of  state  for  a  polytropic  reacting  explosive  mixture 
is  obtained  from  Eq.  (32)  by  setting  »  Kp  *  K. 

We  now  define  the  sound  speed  in  our  reacting  mixture  by  the 
equation 


c 


2 


Kpv 


(33) 


and  we  then  use  Eqs.  (28)  and  (29)  to  express  Eq.  (33)  in  a  more 
convenient  form.  Differentiating  Eq.  (29)  at  constant  s  and  X^  leads  to 
the  equation 


dEs  +  ds  *  0 

x  p 


(34) 


and  it  follows  that  di  =0  because  ds  *  0,  i  »  1,  2,  3  by  assumption. 

P 

Differenting  Eq.  (28)  partially  with  respect  to  p  at  constant  s  and  X.^ 
and  taking  account  of  these  entropy  conditions  gives  the  equation 


11 


ih  ai(1 


(35) 


which  for  our  poly tropic  materials  can  be  rewritten  after  some 
manipulation  as 


K  p 

-P-- 


(v  -  K  Ev  ) 
v  xp  X  ’ 


(36) 


with  K  =  (K  -  K  )/K  .  The  combination  of  Eqs.  (33)  and  (36)  gives 
xp  x  p  x 

the  following  equation  for  K 


K 


5L 


(v  -  K  Ev  ) 
v  xp  X 


(37) 


which  reduces  to  K  =  when  =*  X2  -  X3  -  0  and  v  -  vx>  and  to  K  -  Kp 

when  X^  ■  X^  =*  X3  ■  1.  The  assumption  that  an  explosive  component  is 
released  down  its  own  isentrope  allows  us  to  express  the  Ev^  term  in 
terms  of  pressure  in  a  ZND  wave.  Because  the  explosive  components  are 
polytropic,  their  specific  volumes  can  be  written  as 


(K  -  1) 
o  x 

vx  (K  +  1) 
X 


(38) 


where  6^  *  v“  /vx.  Thus  the  equation  for  £vx  in  our  ZND  wave  can  be 
written  in  terms  of  pressure  as 


Ev 


o  (Kx  _  1)  fpH  //Kx  3  - 
'x  (Kx  +  1)  (  p  ’  iSl  ai(1i  "  V 


(39) 


where  a,  *  a.  v°  /v°  and  ,  E,  a.  *  1. 
i  i  x^  x  i»l  i 
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Nonideal  Steady-State  Detonation  Waves 


The  constitutive  relationship  can  now  be  combined  with  the  Rankine- 
Hugoniot  jump  conditions  and  the  Chapman-Jouguet  conditions  to  generate 
the  equations  governing  the  CJ  state  in  our  nonideal  ZND  wave  and  to 
generate  the  equation  governing  particle  velocity  in  the  steady-state 
reaction  zone. 

Chapman-Jouguet  Conditions 

We  set  (x^ ■  \1  at  the  CJ  point  in  our  steady-state  ZND  wave. 

.AAA 

Then  we  can  represent  ideal  detonation  by  the  set  {X^  =  1,  =  1,  X^  = 

X^  ■  l}  and  nonideal  detonation  by  the  set  { ,  X^,  X^,  X^  }  with  any  one 
of  the  reaction  coordinates  satisfying  the  condition  X  <  1.  In  the 
present  case,  we  assume  that  the  RDX-like  component  X^  and  the  TNT-like 
component  X2  decompose  completely,  and  represent  our  nonideal  detonation 

A  A  A  A 

by  the  set  {X^^  =■  1,  X2  =  1,  X^  <  1,  X^  <  l}.  We  are  therefore  going  to 
consider  incomplete  decomposition  of  the  AN-like  component  and  incom¬ 
plete  reaction  among  the  carbon  from  Xx  and  X2  and  the  oxygen  from  X3. 

We  first  combine  Eqs.  (2),  (4),  and  (33)  to  obtain  the  following 
relationship  between  sound  speed  and  particle  velocity  at  the  CJ  point: 

Ccj  *  K  uCj  (40) 


where 


f1  - 


(41) 


with  Ev  ■  a.(l  -  X,)(v  )  T.  Combining  Eq.  (40)  with  Eqs.  (4)  and  (1) 

X  J  j  w 

then  gives  the  equations 


u 


CJ 


D 

-x - 

K  +  1 


CCJ 


A 

K  D 

A 

K  +  1 
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(42) 


(43) 


K  +  1 


The  equation 


"(3>\ 

U  S'XJ, 


o 

V  —  V 

x  CJ 


expressing  the  tangency  condition  of  the  isentrope  and  the  Rayleigh  line 
at  the  CJ  point  is  then  used  to  obtain  the  following  equations  for  the 
specific  volume  at  the  CJ  point: 


K  +  1  K  +  1 

P  P 


/l  -  K  Sv  /v° ) 
\ _ xp  x  x/ 


K  +  1 
P 


The  combination  of  Eqs.  (44)  and  (47)  then  gives  the  following  equation 
relating  the  polytropic  indices  at  the  CJ  point: 


K  +  1 


K  +  1 


1  -  K  2v  /v 
xp  X  X 


The  equation  relating  the  CJ  pressure  and  the  shock  pressure  at  the 
front  of  our  ZND  wave  is  obtained  as 


2(K  +1)  ^  *  *x/vx) 


by  combining  Eq.  (47),  the  strong  shock  condition 


K  +  1 
x 


and  the  following  equation  expressing  the  balance  of  mass  and  momentum 
in  a  steady-state  wave 

PCJ  (*  '  vCJ/',x) 


(1  -  vH/<) 


The  equation  for  the  detonation  velocity  is  obtained  by  combining  the  RH 
jump  conditions,  the  constitutive  relationship  Eq.  (32),  and  the  CJ 
condition.  It  is  convenient  to  set 

\’i  +  SV>4  ‘  «  <52> 

and  combine  Eqs.  (2),  (3),  and  (32)  to  obtain  the  following  equation  for 
the  particle  velocity: 


u  -  (l  -  K  Ev  /v°)  = 

K  +  1  v  xp  x  x' 

P 


2(K  -  1) 


K  +  1 
P 


The  CJ  condition  D  ■  (K  +  then  gives  the  following  equation  for 

UCJ: 

2  ?(C  _  i \  .  1  2(K  -  1)  A 

UCJ  1  "  K  +  1  ^  “  KxpVTx^  K  +  1  Q 

P  J  P 


where  Q  *  q^  +  q^  +  (q^  +  X^q^)  because  of  our  assumption  that 

X,  ■  X0  *  1.  Combining  Eqs.  (54)  and  (48)  gives  the  following  equation 


1  -  (2K  -  R  1  Ev  /v 

v  xp  xp '  X 


(l  -  K  Ev  /v  ) 
v  xp  x  xJ 


which  with  Eqs.  (44)  and  (48)  leads  to  the  following  equation  for  the 
detonation  velocity: 


2(IT  -  1)Q 
^ _  P 

K  Ev  /v°)  [l  -  (2K 
xp  X  X  L  xp 


K  )Ev  /v°] 

Xp  X  X  J 


The  problem  of  nonideal  detonation  in  our  simple  treatment  of  detonation 

A  A  A 

is  that  of  calculating  E^/v”,  ^3*  ar*d  ^-4*  ft  follows  from  Eq.  (39) 

that  the  equation  for  Ev  /v°  can  be  written  as 

X  X 


E'  v  /v° 

X  X 


(Kx-l)_  .  /PH  \1/Kx 

*  vw  «,(i  -  — 


Kx  +  1  “3(1  -  V  ^pCJ 


but  models  for  calculating  and  X ^  must  be  formulated  or  their  values 
must  be  estimated  before  detonation  parameters  can  be  calculated. 

A  A 

When  a  relationship  between  ant*  ^4  1®  known,  however,  Eqs.  (49), 

(57),  and  (56)  can  be  used  to  investigate  the  dependence  of  D  on  X ^  as 

follows.  We  set  e  *  fl  -  K  Ev  /v°),  write  Eq.  (49)  as 

v  xp  x  x' 

2(K0  +  1)  pCJ 

e  «  - 2 -  (58) 

e  K  +  1  P„  V  3 

x  H 


so  that 


E  v  /v” 

X  X 


-  -  K*  -  1  /"H  \1/K* 

<*3(1  ^3)  ^  4*  X  l  d  / 


1  -  C 


and  assign  e  a  value  a  little  less  than  1.  We  calculate  the  correspond¬ 
ing  values  of  P^j/pfl  and  £vx/vx  from  Eqs.  (58)  and  (59),  and  then  we 
calculate  the  value  of  from  Eq.  (59).  The  relationship  between 

A  A 

and  gives  the  value  of  \^,  and  the  detonation  velocity  can  then 
be  calculated  from  Eq.  (56). 


Examination  of  the  equations  for  the  CJ  parameters  shows  that  the 

A 

equations  for  complete  decomposition  of  X3  with  =  1  have  the  same 
form  as  those  for  polytropic  explosive  with  =  K  =  0.  In  this 
case,  the  equation  for  the  detonation  velocity  reduces  to 


1)2  =  2(Kp  "  +  92  +  q3  +  \q4)  (6°) 

and  the  nonideality  arises  solely  from  the  incompleteness  of  the  reac¬ 
tion  among  the  carbon  from  and  X2  and  the  oxygen  from  Xj- 


Dependence  of  Particle  Velocity  on  Reaction  Coordinates 

It  is  convenient  to  use  Eq.  (53)  to  derive  the  relationship  of  the 
particle  velocity  and  the  reaction  coordinates  in  a  steady-state  ZND 
wave.  We  write  Eq.  (53)  to  make  its  left-hand  side  into  a  perfect 
square  and  then  combine  this  equation  with  the  CJ  conditions,  Eqs.  (42) 
and  (48),  to  obtain  the  equation 


(l  -  K  Ev  /v°) 
u  v  xp  X  x' 

1  + 

,  .  2(KP  -  11Q 

1/2- 

UCJ  (l  -  K  Ev  /v°) 

v  xp  X  X  ' 

<l  -  Wv°)v 

Equation  (56)  then  allows  Eq.  (61)  to  be  written  as 


CJ 


(l  -  K  Ev  /v° ) 

(l  -  K  Ev  /v° ) 
v  xp  X  X ' 


1  + 


1  - 11  -  (2Kxp  -  ^xpK^1  - 


o  \2 


(l  -  K  Ev  /v  ) 
v  xp  X  X  ’ 


l/2> 


/( 62) 
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* 


iSl  at(l  -  V 


(63) 


Sv  /v 

X  X 


/UH  ^ 

f/KX  1 

^ucj\ 

l^cjy 

from  Eq.  (39),  and  with  =  1,  =  X^,  and  X^  =  X^  at  the  CJ 

point  where  Q  =  Q,  Evx  =  Evx>  and  u  =  uqj.  The  equation  for  the 
particle  velocity  in  a  steady-state  ZND  wave  when  X^  decomposes 

A, 

completely  is  obtained  by  setting  X^  *  1  in  Eq.  (62)  as 


-S—  -  (1  -  K  Zv  h°)  1  +  1  -  - 2 - — 3;—*— 

UCJ  *P  X  X  (1  -  K  Sv  /V°)2  Q  I 

xp  x  x '  d 


where  Qd  =  9^  +  92  +  q3  +  ^4q4*  7116  e9uation  f°r  the  Particle 

velocity  in  a  steady-state  ZND  wave  in  polytropic  explosive  is  obtained 
from  Eq.  (62)  by  setting  *  K  =  0  and  is  conveniently  written  as 

r  „  „  ~  « ii/2 

=  1  +^(1  -  \)  +  a2d  “  V  +  V1  '  W  +  a4(1  "  W  (65) 


where 

* 

aL  =  q  x /Q 
«2  ■  q2/Q 
“3  =  X3q3/2 

®4  =  X494/Q 


al  +  a2  +  a3  +  a4  3  * 


Differentiating  Eq.  (53)  with  respect  to  time  along  a  particle  path 
gives  the  following  equation  for  the  deceleration  in  a  steady-state 
reaction  zone: 


Because  3u/3t  is  finite,  it  follows  in  general  from  Eqs.  (66)  and  (67) 
that  the  rates  of  all  the  reactions  must  be  zero,  i.e.,  ax^/dt  =  0  for  i  = 

1,  2,  3,  and  4,  at  the  sonic  point  where  u  =  uqj  =  o(l  ~  ^xp^x^x^ (Kp  +  0* 
In  other  words,  the  CJ  point  is  located  at  the  end  of  the  steady-state 
reaction  zone  of  our  nonideal  detonation  wave.  The  following  equation, 


UH  \at7 


2(K  -  1)<3U  K  (K  -  1)  «  ,  ax, 

p  H  ^  xp'  x  '2  -  i 

K  +  1  +  k  +  i  ^  i«l  “i  atT 

P  P  H 


obtained  by  setting  u  “  ujj  and  X^  =  X2  =  X^  =  X^  =  0  in  Eq.  (66),  leads 

to  the  conclusion  that  the  particle  velocity  will  Increase  at  the  front 

of  the  wave  when  K  (K  -  1)  a,  3X, /3tu  >  2(K  -  1)  (cJu/u^).  Waves 

xp  x  i*l  i  i  H  p  H' 

with  this  property  will  not,  however,  be  considered  in  this  report. 
Equations  (62),  (64),  and  (65)  can  be  used  to  calculate  Lagrange  par¬ 
ticle  velocity  histories  when  expressions  for  the  reaction  coordinates 
are  known  as  functions  of  time.  The  Lagrange  histories  can  then  be 
examined  to  investigate  the  influence  of  the  relative  rates  of  the 
exothermic  reactions  on  particle  velocity.  Such  calculations  will  be 
presented  in  the  next  section. 

♦ 
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Model  Nonideal  Detonation  Calculations 

Properties  of  the  explosive  X,  its  components  X^>  x2,  and  X3,  and 
their  detonation  products  were  chosen  to  model  Amatex  20.  We  set 
v°  m  0.6201  cm^/g  so  that  the  initial  density  p*  ■  1.613  g/cm^,  and  we 
assigned  X^,  X2,  and  X3  the  values  of  the  parameters  shown  in  Table  1. 

A 

We  set  *  587.5  cal/g  and  q^  -  235  cal/g  according  to  Eqs.  (23) 
and  (24),  and  we  set  ^  =  5.34  and  Kp  *  2.67,  so  that  *  0.6152  and 
*  0.5.  In  this  case,  q^  =  i£i»  9^  ”  951.1  cal/g,  q^  +  q^  ■  1186.1 
cal/g,  and  when  Xlf  X2,  and  X3  decompose  completely,  the  detonation 
velocity  D,  the  shock  pressure  p^,  the  CJ  pressure  p^j,  and  the  CJ 
particle  velocity  uCj  vary  from  values  of  6.98  mm/ps,  247.7  kbar,  214 
kbar,  and  1.90  mm/ps  to  values  of  7.80  mm/ps,  309.6  kbar,  267.4  kbar, 
and  2.12  mm/ps  as  A ,  varies  from  zero  to  one.  We  will  now  calculate 

^  A  A  A  A 

detonation  parameters  when  A^  »  X^  *  1,  X^  <  1,  and  X^  <  1. 

Detonation  Parameters 

We  calculate  detonation  parameters  for  a  nonideal  detonation  with 

A  A 

X^  >  1  and  X^  <  1.  Our  calculations  are  based  on  the  intuitive  hypo- 

A  a 

thesis  that  the  slope  of  the  D  versus  A3  curve  is  positive  as  X^  approa 
ches  1  when  X^  is  a  function  of  We  accordingly  differentiate  Eq. 

(56)  and  set  dD/dX^  >  0  when  X3  -  1  to  find  conditions  to  satisfy  this 
hypothesis.  For  notational  simplicity,  we  set  R  *  and  21?^  -  ■ 

K2*  Differentiating  Eq.  (56)  with  respect  to  A^  then  gives  the  equation 


2  dD 


1  dQ 

T  - JT"  T 

Q  <1X3 


which  reduces  to 


2  dD 

- x — 

D  dX3 


d(S\/vx) 

h  +  2lti 

dX3 

(1  -  -  W\). 

q3  +  q4^X4  +  dX4/dX3) 

ql  +  q2 

A 

+  <J3  +  VU 

A 

A~  ■  1 

(K  -  1)  /2(K  +  l)\1/Kx 


(69) 


when  X^  “  1*  We  set  X^  =  n  X^  with  n  <  1  for  simplicity  so  that 
(dD/dX,)^  >  0  when 

X3=l 

q3  +  2n  q  _  (K  -  1)  /2(K  +  l)\1/Kx 

qL  +  q2  +  q3  +  n  ^4  >  2  Kxp  013  (Kx  +  1)  \  (KX  +  1) j  (71) 

Solving  this  inequality  for  n  shows  that 
(dD/dXo  )  A  >  0 

X3  "  1 

when  n  >  0.454.  We  are  now  in  a  position  to  calculate  detonation 

A 

parameters  for  specified  values  of  n  and  X^-  We  set  n  =  0.6,  and  use 
Eqs.  (49),  (57),  (56),  (48),  (42),  (2),  and  (46)  to  calculate  values  of 

A 

the  detonation  parameters  for  specified  values  of  X~  as  follows. 

J  A 

Equations  (49)  and  (57)  are  solved  for  a  given  value  of  X~  to  find  the 

A  ^ 

corresponding  values  of  (py/pcJ^  anc*  Ev„/V  •  Equations  (56)  and  (48) 
are  then  used  to  calculate  the  values  of  D  and  K,  and  equations  (42), 

(2),  and  (46)  are  used  to  calculate  the  values  of  u^j,  pqj  and  vqj/v°. 
Calculated  values  of  the  detonation  parameters  are  given  in  Table  2. 

Values  of  the  detonation  velocity  are  given  to  three  figures  in 
Table  2  to  show  that  the  detonation  velocity  has  a  minimum  at  X^  *»  0.75, 
while  the  CJ  pressure  and  particle  velocity  increase  monotonically  as  X^ 
increases  from  a  value  of  0.45  to  a  value  of  1.  The  detonation  velocity 

A 

exhibits  a  minimum  because  K  decreases  monotonically  with  X^  as  Uqj 
increases.  This  minimum  demonstrates  the  influence  of  the  equations  of 
state  of  the  condensed  explosives  on  the  nonideal  detonation  process  in 
composite  explosives. 
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Lagrange  Particle  Velocity  Histories 


It  is  clear  from  Eq.  (62)  that  expressions  for  time  variations  of 


the  reaction  coordinates  ,  X„ ,  X, ,  and  X.  must  be  known  before  we  can 

1  2  3  4 

calculate  the  Lagrange  particle  velocity  history  in  the  reaction  zone  of 
a  nonideal  steady-state  reaction  zone.  Bearing  in  mind  that  the  rates 
of  all  the  exothermic  reactions  nust  be  zero  at  the  CJ  point,  we  formu¬ 
lated  simple  expressions  for  the  reaction  coordinates  to  perform  such 
calculations.  These  expressions  are  based  on  the  following  equation  for 
a  reaction  that  proceeds  to  completion. 

-in. 


(i  -  xt)  • 


1  - 


(t  -  -c) 


(72) 


where  n^  >  2,  the  Lagrange  time  x  denotes  the  time  a  particle  enters  the 
wave,  and  the  overall  reaction  time  T  denotes  the  time  it  takes  a 
particle  to  travel  from  the  shock  front  to  the  CJ  plane.  Equation  (72) 
is  written  so  that  X^  =  0  at  the  shock  front  where  t  -  x  =  0,  and  X.^  =  1 
at  the  CJ  point  where  (t  -  x)  =  T.  Differentiating  Eq.  (72)  with 
respect  to  time  along  a  particle  path  gives  the  equation  for  the 
reaction  rate  as 


and  it  follows  from  Eq.  (73)  that  we  are  restricting  our  consideration 
to  reactions  with  negligible  activation  energies. 


We  used  Eq.  (72)  with  i  =  1  and  i  =  2  to  model  respectively  the 
decomposition  of  Xj^  and  %2’  To  account  for  the  fact  that  the  faster 
reacting  RDX-like  component  X^  may  decompose  completely  within  the 
reaction  zone,  we  introduce  a  parameter  f  >  1  and  modify  Eq.  (72)  by 
setting 


1 


1 

m 


0  <  (t  -  t)/T  <  (f  -  l)/f 


(74) 
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and 


X  =  1  for  (t  -  t)/T  >  (f  -  l)/f  (75) 

In  this  case,  the  decomposition  of  X^  is  complete  when  (t  -  x)/T  * 

(f  -  l)/f;  when  f  =  3,  for  example,  the  Xj^  decomposition  time  is  2T/3. 

It  is  clear  that  Eq.  (74)  can  be  written  in  the  same  form  as  Eq.  (72)  by 
scaling  the  reaction  time  and  setting  the  new  reaction  time 
as  T  -  (f  -  1)  T/f. 

To  model  the  case  when  the  decomposition  of  X^  is  not  self- 
propagating,  we  simply  set 


3X3  9X2 

W~  m  n3  at- 


(76) 


with  n^  <  1  and  ax^/at  given  by  Eq.  (73)  with  i  *  2,  so  that  the 
decomposition  of  X^  ceases  at  the  end  of  the  reaction  zone  where  it  is 
not  being  supported  by  the  decomposition  of  the  other  more  energetic 
components.  Integrating  Eq.  (76)  gives  the  equation 


X3  =  n3  X2  (77) 

and  it  follows  that  X^  will  have  a  value  of  *  n^  at  the  end  of  the 

A 

reaction  zone  where  =  1. 

To  model  the  case  when  the  recombination  reaction  of  the  oxygen 
from  X3  and  the  carbon  from  X!  and  X2  proceeds  at  a  lower  rate  than  the 
production  of  the  oxygen,  we  set 

ax4  axq  ax2 

nr  "  n4  ”5t~  n4n3  at-  (78) 


with  n^  <  1,  so  that 

x4  3  n^X^  a  n4n3X2 


(79) 
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and  X^  *  n^X^  *  n4n3  at  t"e  en<*  t*ie  reacti°n  zone  where  *  1. 
Equations  (72),  (74),  (75),  (77),  and  (79),  with  Eq.  (53)  rewritten  as 


n^n^  at  the  end  of  the  reaction  zone  where 


— V  ■  Af— ' 

UCJ  J  ^UCJ 


fl  -  K  Ev  /v°l 
1  xp  X  x ; 


1  -  A  ^  -  Kxp&x/vxl  | 


with  A  *  2D/(Kp  +  1)  u^j,  were  used  to  calculate  Lagrange  particle 
velocity  histories  in  the  reaction  zone  of  nonideal  detonation  waves 

A  A 

with  X ^  =  0.8  and  X^  *  0.6.  In  this  case,  D  =  7.46  mm/ps,  Uqj  =  1.977 


mm/ps,  and  u  /u  T  ■  1.190;  the  equations  for  calculating  values  of 
H  CJ 

X3  and  X4  are  obtained  as  X3  =*  0.8  X2  and  X4  -  0.6  X3  by  setting  n3  = 

0.8  and  n,  =  0.6  in  Eqs.  (77)  and  (79). 

4 

A  Newton-Raphson  method  was  used  to  calculate  values  of  (u/u^j) 
from  Eq.  (88)  throughout  the  reaction  zone  with  values  of  X^  and  X^ 
calculated  from  Eqs.  (72),  (74),  and  (75).  The  results  of  three  such 
calculations,  based  on  the  assumption  that  X^  decomposes  faster  than  X£, 
are  shown  graphically  in  Figures  la-3a  and  lb-3b.  The  profiles  in 
Figures  la  and  lb,  2a  and  2b,  and  3a  and  3b  were  calculated  respectively 
with  the  following  sets  of  parameters:  {n^  =3,  f  •  4,  *  2}, 

fnx  -  2,  f  =  3,  n2  »  2},  and  =  3,  f  -  3,  n2  -  3}. 

Differentiating  Eq.  (66)  partially  again  with  respect  to  t  shows  in 
general  that  (Qu/9t)CJ  <  0  when  f  >  0,  and  that  (du/dt)Cj  m  0 

when  ( d 2  X^dt2)^  ~  0  when  X1 ,  X3  and  X^  are  governed  by  Eqs.  (74),  (76), 
and  (78).  Examination  of  the  particle  velocity  histories  shows  that  three 
conditions  are  satisfied  because  (&u/5t)rT  <  0  in  Figures  la  and  2a  where 
n2  ■  2  and  ( &  X^/ >  0,  and  (du/ftt)^  =  0  in  Figure  3a  where  n£  51  3 
and  (d^X2/3t^)Cj  *  0.  Further  examination  of  Figures  la-3a  show  that 
detailed  properties  of  the  Lagrange  particle  velocity  history  depend  on 
the  equations  of  state  of  the  explosive  and  the  detonation  products  as 
well  as  the  energy  release  rate.  However,  the  qualitative  features  of 
the  particle  velocity  history  are  governed  by  the  energy  release  rate, 
as  shown  by  the  observations  that  (1)  the  history  calculated  with  (n^  *  3, 
no  ■  2)  is  generally  steeper  than  the  history  calculated  with  (nj  =  2, 
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FIGURE  la  PARTICLE  VELOCITY  vs  TIME  IN  REACTION  ZONE  CALCULATED 
WITH  n,  =  3,  f  =  4,  AND  n2  =  2 
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DIMENSIONLESS  PARTICLE  VELOCITY,  u/u 
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FIGURE  1b  REACTION  COORDINATES  vs  PARTICLE  VELOCITY  IN  REACTION  ZONE 
CALCULATED  WITH  n,  =  3,  f  =  4,  AND  n,  =  2 


DIMENSIONLESS  TIME,  It-TAUJ/T 


JA-1 290-9 


FIGURE  2a  PARTICLE  VELOCITY  vs  TIME  IN  REACTION  ZONE  CALCULATED 
WITH  n,  =  2,  f  =  3,  and  n-  =  2 
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FIGURE  2b 


CJ 


REACTION  COORDINATES  vs  PARTICLE  VELOCITY  IN  REACTION 
ZONE  CALCULATED  WITH  ru  =  2,  f  =  3,  and  n„  =  2 


DIMENSIONLESS  PARTICLE  VELOCITY,  u/u. 
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FIGURE  3b 
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REACTION  COORDINATES  vs  PARTICLE  VELOCITY  IN  REACTION  ZONE 
CALCULATED  WITH  n,  =  3,  f  =  3.  AND  n„  =  3 


n2  =  2),  and  (2)  the  history  calculated  with  (np  =  3,  n2  =  2)  is 
generally  steeper  than  the  history  calculated  with  (np  =3,  r\2  ~  2). 

Conclusions 

We  have  formulated  a  treatment  of  steady-state  detonation  waves  in 
composite  explosives.  Our  formulation  extends  the  classical  ZND  model 
for  ideal  detonation  into  a  prototype  model  for  nonideal  detonation. 

The  consititutive  relationship  used  to  describe  shocked  reacting  explo¬ 
sive  is  based  on  the  assumption  that  the  explosive  and  its  detonation 
products  are  governed  by  different  mechanical  equations  of  state  and 
attain  mechanical,  but  not  thermal,  equilibrium.  The  kinetic  processes 
responsible  for  nonideal  behavior  are  assumed  to  be  incomplete  decom¬ 
position  of  the  slowest  reacting  explosive  component  and  the  incomplete 
attainment  of  chemical  equilibrium  among  the  detonation  products  from 
the  different  explosive  components. 

To  make  our  model  applicable  to  Amatex  20,  we  considered  a  compo¬ 
site  explosive  containing  an  RDX-like  component  Xp>  a  TNT-like  component 
x2>  and  an  AN-like  component  X3.  In  this  case  the  constitutive 
relationship  must  include  four  reaction  coordinates:  X^,  X^,  and  X^ 
describe  the  decomposition  of  Xp,  X2,  and  X3;  X4  describes  the  reaction 
of  the  oxygen  produced  by  X3  and  the  carbon  produced  by  Xp  and  X2*  The 
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reaction  coordinates  at  the  end  of  the  reaction  zone,  X^,  \^>  X3>  and  X4, 
satisfy  the  conditions  X^  =  X2  =  X^  =  X4  =  1  in  an  ideal  detonation 

AAA  a 

wave,  and  satisfy  the  conditions  X^  =  X2  =  1 ,  X^  <  1,  and  ^  <  1  in  a 
nonideal  detonation  wave. 

The  consititutive  relationship  was  combined  with  the  Rankine-Hugoniot 
jump  conditions  and  the  CJ  sonic  condition  to  obtain  the  equations  for 
the  CJ  parameters  at  the  end  of  the  reaction  zone  in  a  steady-state 
nonideal  detonation  wave.  The  equation  governing  particle  velocity  was 
combined  with  the  CJ  condition,  relating  particle  velocity  and  detona¬ 
tion  velocity  in  a  nonideal  detonation  wave,  to  obtain  the  equation 
relating  the  detonation  velocity  to  the  chemical  energy  released  in  the 
wave.  These  equations  for  nonideal  detonation  specify  the  dependence  of 


the  detonation  parameters  on  the  reaction  coordinates  and  demonstrate 
that  the  detonation  process  is  influenced  by  the  equation  of  state  of 
the  explosive.  Detonation  calculations,  using  explosive  parameters 
chosen  to  model  Amatex  20,  were  performed  with  values  of  X^  the  range 

a  a  ” 

0.45  <  X^  <1  and  X^  =  0.6  X^  to  investigate  the  dependence  of  nonideal 
detonation  parameters  on  the  degree  of  incomplete  decomposition  of  the 
explosive.  Simplistic  reaction  rate  functions  were  also  constructed  and 
used  to  calculate  Lagrange  particle  velocity  histories  in  nonideal 
detonation  waves. 

We  conclude  that  the  result  of  our  theoretical  study  is  a  well- 
formulated  model  for  steady-state  detonation  that  provides  us  with  a 
better  understanding  of  nonideal  detonation.  However,  we  are  still 
faced  with  the  major  problem  of  formulating  a  realistic  treatment  of  the 
decomposition  of  AN  and  subsequent  reactions  among  the  detonation 
products.  Before  the  present  model  of  nonideal  detonation  can  be  used 
to  calculate  realistic  nonideal  detonation  parameters,  we  must  incor¬ 
porate  (1)  models  for  these  kinetics  processes  that  predict  values  of 

A  A 

X-  and  X,;  and  (2)  more  realistic  equations  of  state. 


Ill  EXPERIMENTS 


Basic  Technique 

The  Lagrange  gage  and  analysis  method  used  In  this  work  is  based 
on  the  experimental  measurement  of  Lagrangian  particle  velocity 
histories  within  an  explosive  target  undergoing  planar  detonation.  Our 
experimental  techniques  have  been  described  previously*’ * **■ » ^  and  are 
summarized  here. 

The  experimental  configuration  is  shown  in  Figure  4.  A  target  in 
the  form  of  a  right  circular  cylinder,  typically  about  0.2  m  (8  inches) 
in  diameter  and  0.13  m  (5  inches)  high,  is  constructed  from  the  explo¬ 
sive  to  be  characterized.  A  plane-wave-initiated  high  explosive  charge 
loads  the  target  with  a  pulse  of  the  desired  amplitude  and  shape. 
Particle  velocity  histories  are  recorded  at  several  Lagrange  positions 
(depths)  in  the  target  for  several  microseconds  using  embedded  electro¬ 
magnetic  particle  velocity  (EPV)  gages.  Expendable  Helmholtz  coils 

wound  on  plywood  frames  provide  a  uniform  quasi-static  uniaxial  magnetic 
induction  in  the  region  of  the  EPV  gages.  The  gages,  targets,  and  data 
are  described  below. 

EPV  Gages 

Basic  Operation 

EPV  gages  are  based  on  the  principle  that  an  electrical  conductor 
moving  in  a  fixed  magnetic  field  will  generate  a  motional  electromotive 
force  (EMF )  proportional  to  its  velocity.  If  the  conductor,  motion,  and 
magnetic  induction  are  each  uniaxial  and  mutually  orthogonal,  as  in  the 
present  case,  then  by  Faraday's  law  of  induction  for  moving  circuits 


The  term  Lagrange  refers  to  an  observation  made  at  a  specific  material 
element  in  the  flow. 
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FIGURE  4  BASIC  CONFIGURATION  OF  LARGE-SCALE  LAGRANGE 
PARTICLE  VELOCITY  GAGE  EXPERIMENTS 

The  vertical  lines  in  the  target  represent  cross  sections  of  the 
active  elements  of  individual  embedded  electromagnetic 
particle  velocity  gages;  ten  gages  are  indicated. 


where  E  is  the  EMF,  t  is  time,  B  is  the  magnitude  of  the  magnetic  induc¬ 
tion,  H  is  the  length  of  the  conductor,  and  u  is  the  velocity  of  the 
conductor. 

The  EPV  gages  used  in  this  work  were  formed  from  aluminum  strips 
nominally  0.15  mm  thick  and  3  mm  wide.  The  strips  were  bent  into  the 
rectilinear  U-shape  shown  in  Figure  5.  The  crossbar  of  the  U  is  the 
active  element  of  the  gage  and  is  nominally  25.4  mm  long.  The  sides  of 
the  U  serve  as  electrical  leads  that  carry  the  signal  out  the  back  of 
the  target;  these  leads  do  not  contribute  to  the  signal  because  they  are 
strictly  parallel  to  the  direction  of  motion  in  the  target.  As  in  our 
previous  work,^»^^»^^  the  gages  were  not  insulated  because  their  cross- 
sectional  area  was  assumed  to  be  great  enough  to  result  in  a  low 
internal  gage  resistance  compared  with  alternative  conductive  paths  in 
the  reaction  products.  After  a  ringup  time  of  about  0.1  ps,  the  gages 
are  estimated  to  reach  mechanical  equilibrium  with  the  surrounding  flow 
so  that  u(t)  is  equal  to  the  local  particle  velocity. 

The  gage  outputs,  E(t),  were  determined  from  the  potential  drop 
across  viewing  resistors  connected  to  the  gages  in  simple  series 
circuits.  Standard  high-frequency  circuit  components,  consisting  of 
shielded  RG213  or  RG58  coaxial  cables,  the  viewing  resistors,  and 
Tektronix  7000  Series  oscilloscopes  or  Nicolet  digital  oscilloscopes, 
were  used  in  measuring  the  potential  drops.  The  oscilloscopes  were 
operated  in  normal  (single-ended),  rather  than  differential,  mode. 
Because  the  resistance  of  the  gages  and  cables  was  small  compared  with 
the  viewing  resistance,  E(t)  was  taken  to  be  equal  to  the  recorded 
potential  drops.  In  this  work  B  was  about  450  gauss,  and  the  peak 
particle  velocities  were  approximately  2  mm/ps,  so  E[nax  was  greater  than 
2  volts  and  presented  no  recording  problems.  Analog  oscilloscope  time- 
bases  were  operated  at  either  500  or  1000  ns/cm  and  were  time-correlated 
to  less  than  50  ns.  Oscilloscope  voltage  and  time  calibrations  were 
checked  at  the  beginning  and  conclusion  of  each  shot  sequence. 
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FIGURE  5  GAGE  BLOCK  USED  TO  CAST  EPV  GAGES  INTO 
EXPLOSIVE  TARGETS 

The  gage  block  shown  contains  ten  gages. 

The  configuration  of  the  target  to  be  cast  around  it 
is  indicated  by  the  dashed  lines. 


Emplacement 

The  gages  were  mounted  on  blocks  and  embedded  in  the  targets  during 
casting,  as  indicated  in  Figure  5.  The  gage  blocks,  machined  from 
linen-phenolic,  each  contained  ten  gages.  The  blocks  were  secured  to 
the  bottom  of  cylindrical  molds,  and  the  explosive  was  cast  around  the 
assemblies  so  that  the  gages  were  not  disturbed.  Gage  positions  in  the 
completed  targets  were  checked  radiographically,  and  no  displacements 
were  detected.  Gage  stability  during  explosive  casting  was  also  demon¬ 
strated  in  a  previous  program  in  which  a  gage  block  was  embedded  in  an 
explosive  simulant  mixed  to  the  efflux  viscosity  of  an  Amatex  20  explo¬ 
sive  and  then  destructively  sampled. ^  Casting,  machining,  and  x-raying 
of  the  explosives  were  performed  in  the  M-l  Group  of  Los  Alamos  National 
Scientific  Laboratories  (LANSL)  under  the  direction  of  A.  P.  Torres  and 
J.  Sanchez. 


Targets 

Coarse-Grained  Amatex  20 

The  variation  of  composition  and  density  in  the  coarse-grained 
Amatex  20  targets  was  determined  at  LANSL  from  corings  of  a  quality- 
control  casting. 

As  reported  by  LANSL,  the  coring  pattern  is  shown  in  Figure  6,  and 
the  Amatex  20  composition  and  deroity  and  the  ammonium  nitrate  prill 
particle  size  are  given  in  Tables  3  and  4. 

Target  Configuration 

Standardized  gage  block  assemblies  and  target  dimensions  were  used 
in  the  Lagrange  gage  experiments.  Figure  7  shows  the  configuration  and 
the  gage  numbering  system  and  gives  the  nominal  target  size  and  gage 
block  position  to  the  nearest  millimeter. 

The  gage  locations  in  the  axial  (wave  propagation)  direction  were 
controlled  to  high  precision.  Each  block  contained  two  sets  of  five 
identically  positioned  gages  (numbered  1  through  5  and  6  through  10  as 


1 

19.47 

41.22 

39.31 

1.607 

2 

21.19 

42.04 

36.77 

1.571 

3 

18.67 

41.89 

39.44 

1.603 

4 

17.62 

42.40 

39.98 

1.572 

Average: 

19.24 

41.88 

38.88 

1.588 

Spread: 

3.57 

1.18 

3.21 

0.036 

Table  4 

AMMONIUM  NITRATE  PRILL  PARTICLE  SIZE, 
COARSE-GRAINED  AMATEX  20 


Mesh  Size 

(m®) 

Cumulative  Percent 
Retained 

2800 

0 

2000 

22.0 

1700 

62.4 

1400 

88.3 

1000 

98.7 

710 

100.0 
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shown  in  Figure  7).  Thus  there  were  two  gages  in  each  of  five  planes 
parallel  to  the  initiating  plane.  The  nominal  distance  from  the 
initiating  plane  (front  surface  of  target)  to  the  first  gage  plane  was 
24  mm  (0.93  inch);  the  distance  between  successive  gage  planes,  used  in 
the  analysis,  was  12.7  ±  0.1  mm.  The  lateral  distance  (center  to 
center)  between  adjacent  gages  was  12.0  mm  ±-  0.2  mm. 

Data 

Two  experiments  were  performed  with  coarse-grained  Amatex  20.  A 
203-mm-diameter  (8-inch)  P-80  plane  wave  lens  in  contact  with  a  1-inch 
by  8-inch  pad  of  Composition  B  was  used  to  initiate  the  Amatex  20  target 
in  the  first  experiment.  A  203-mm-diameter  (8-inch)  P-80  plane  wave 
lens  in  contact  with  a  1-inch  by  8-inch  pad  of  TNT  was  used  to  initiate 
the  Amatex  20  target  in  the  second  experiment.  Particle  velocity 
records  were  obtained  in  both  experiments.  All  high-resolution  gage 
records  except  5  and  6  were  obtained  with  Tektronix  7844  dual-beam 
oscilloscopes;  high  resolution  records  for  gages  5  and  6  were  recorded 
on  Tektronix  7704A  single-beam  oscilloscopes.  The  oscilloscope  sensi¬ 
tivities  for  the  high-resolution  records  were  500  mV/cm,  vertical,  and 
500  ns/cm,  horizontal.  A  beam  intensity  modulation  system  was  used  to 
simultaneously  generate  several  short  blanks  on  all  the  high  resolution 
oscilloscope  records.  The  blanks  permited  the  individual  gage  records 
to  be  correlated  on  an  absolute  time-scale. 

The  gage  records  obtained  in  the  first  experiment  with  the 
Composition  B  driver  are  not  included  here  because  they  show  that  the 
detonation  wave  in  coarse-grained  Amatex  20  was  overdriven.  The  as- 
recorded  EPV  gage  profiles  obtained  in  the  second  experiment  with  the 
TNT  driven  are  shown  in  Figure  8.  These  records  are  the  basis  for  the 
Lagrange  analysis  presented  in  the  next  section. 

Two  features  of  the  records  pertaining  to  the  experimental 
procedures  are  discussed  here.  First,  high  frequency  oscillations  are 
present  on  most  gage  records  during  or  just  after  the  initial  accelera¬ 
tion  of  that  gage.  The  amplitude  of  the  oscillations  decreases  with 
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depth  in  the  target  from  the  initiating  plane.  Because  the  cause  is 
unknown  and  the  oscillations  do  not  interfere  with  data  analysis  except 
in  the  first  two  gage  planes,  the  phenomenon  was  simply  labeled  noise 
and  was  not  studied.  Factors  contributing  to  this  noise  may  be  the 
length  of  the  gage  leads  (which  may  serve  as  antennas)  and  the  conduc¬ 
tivity  of  the  explosive,  especially  in  the  reaction  zone;  the  phenomenon 
is  probably  not  caused  by  a  real  flow  oscillation. 

The  second  feature  to  note  is  the  disturbance  on  a  particular  gage 
record  when  downstream  gage  enter  the  wave  front.  Again,  we  believe 
that  these  disturbances  are  electrical  in  origin  and  not  the  result  of 
real  flow  oscillations.  The  disturbances  are  not  large  enough  to  inter¬ 
fere  with  data  analysis,  and  they  are  useful  for  time-correlating  the 
records.  They  may  be  caused  by  resistive  coupling  of  gage  circuits 
through  the  highly  conductive  wave  front/reaction  zone. 


Lagrange  Analysis 
Procedure 

Lagrange  analysis  based  on  particle  velocity  provides  a  quantitative 
description  of  the  one-dimensional  flow  produced  in  a  shock  or  detonation 
wave  experiment  from  a  set  of  particle  velocity  gage  histories  recorded 
along  different  particle  paths.  The  flow  variables  are  generated  by 
integrating  the  equations  of  motion  with  flow  derivatives  estimated  from 
the  data  supplied  by  the  gage  records.  The  momentum  equation 


wr  ■  fc-  <8 

h  t 

where  t  and  h  denote  time  and  Lagrange  distance,  and  the  subscript  o 
denotes  the  unshocked  conditon,  gives  dynamic  Lagrange  (p,u)  relation¬ 
ships  between  pressure  and  particle  velocity.  The  continuity  equation 


dv  .to 


(83) 
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gives  dynamic  Lagrange  (v,u)  relationships  between  specific  volume  and 
particle  velocity.  The  elimination  of  v  between  these  (p,u)  and  (v,u) 
relationships  gives  the  corresponding  Lagrange  (p,v)  relationships. 


The  voltage-time  oscilloscope  traces  recorded  by  the  particle 
velocity  gages  in  detonating  coarse-grained  Amatex  20  (Figure  8)  were 
the  basis  for  the  Lagrange  analysis  performed  here.  Before  the  experi¬ 
mental  data  can  be  used  in  a  Lagrange  analysis,  the  oscilloscope  traces 
must  be  digitized  and  converted  into  particle  velocity  histories.  The 
particle  velocity  histories  were  given  the  same  numbers  as  the  gages 
from  which  they  were  generated.  It  follows  from  Figure  7  that  particle 
velocity  histories  were  recorded  by  pairs  of  gages  in  each  of  five 
Lagrange  planes  in  the  sequence  (1,10),  (2,9),  (3,8),  (4,7),  and  (5,6). 
Selection  and  preparation  of  a  single  set  of  particle  velocity  histories 
for  Lagrange  analysis  is  discussed  next. 

Selection  and  Smoothing  of  Particle  ;locity  Records 

First,  the  present  particle  velocity  histories  recorded  in 
detonating  Amatex  20  containing  coarse-grained  AN  were  compared  with 
those  recorded  earlier^1)  in  detonating  Amatex  20  containing  baseline 
AN.  The  histories  recorded  in  the  coarse-grained  Amatex  20  were  found 
to  be  more  noisy  than  those  recorded  in  the  baseline  Amatex  20  as 
expected.  The  records  produced  by  gage  pairs  (1,10)  and  (2,9)  were 
considered  to  be  too  noisy  for  analysis  and  were  discarded.  The 
oscilloscope  traces  recorded  by  the  pairs  of  gages  (3,8),  (4,7),  and 
(5,6)  were  digitized  and  converted  into  particle  velocity  histories. 

We  first  compared  the  particle  velocity  histories  recorded  by  pairs 
of  gages  located  in  the  same  Lagrange  plane.  This  comparison  showed 
that  the  set  of  particle  velocity  histories  exhibit  the  following 
characteristics : 

•  Histories  3  and  8  do  not  agree  for  approximately  0.5  ps 
at  the  wave  front  but  are  in  good  agreement  over  the 
last  3  pa. 


•  Histories  4  and  7  are  noisy  for  »  0.2  ps  at  the  wave 
front  but  are  in  excellent  agreement  over  ~  4  ps. 

•  Histories  5  and  6  are  in  very  good  agreement  over 
~  4  J13 

We  then  compared  the  particle  velocity  histories  (3,  4,  5,)  from  one 
side  of  the  charge,  then  the  corresponding  particle  velocity  histories 
(8,  7,  6)  from  the  other  side  of  the  charge,  and  finally  all  the 
particle  velocity  histories  against  each  other. 

The  set  (3,4,5)  showed  that 

•  Histories  3,  4,  and  5  have  respective  peak  particle 
velocities  of  ~  2.2  mm/ps,  ~  2.3  mm/ps  and  2.0  mm/ ps  at 
the  wave  front. 

•  Apart  from  exhibiting  different  peak  structures, 
histories  4  and  5  agree  for  ~  0.6  ps  after  the  arrival 
of  the  wave  front. 

•  The  relative  positions  of  histories  3,  4,  and  5  over  the 
last  ~  3.4  ps  are  compatible  with  the  flow  in  a  Taylor 
wave. 

The  set  (8,  7,  6)  showed  that 

•  Histories  8,  7,  and  6  have  respective  peak  particle 
velocities  of  ~  2.3  mm/ps,  ~  2.3  mm/ps,  and  2.0  mm/ps. 

•  Apart  from  history  6  exhibiting  a  different  peak 
structure  than  histories  8  and  7,  the  agreement  among 
these  histories  for  ~  0.6  ps  indicates  the  presence  of  a 
steady-state  reaction  zone  in  coarse-grained  Amatex  20. 

•  The  relative  positions  of  histories  8,  7,  and  6  over 
their  last  ~  3.4  ps  are  compatible  with  the  flow  in  a 
Taylor  wave. 

Comparison  of  all  the  particle  velocity  histories  showed  that 

•  Histories  4,  8,  and  7  are  in  good  agreement  after  the 
arrival  of  the  wave  front  for  ~  0.6  ps  as  the  particle 
velocity  drops  to  ~  1.6  mm  ps,  and  except  for  a 
difference  in  peak  structure  the  histories  also  agree 
with  histories  5  and  6  over  the  same  time  interval. 

The  relative  positions  of  the  records  over  their  last 
~  3.4  ps  are  compatible  with  the  flow  in  a  Taylor  wave. 
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We  concluded  that  there  are  two  distinct  regions  of  flow  behind  the  wave 
front  in  detonating  coarse-grained  Amatex  20:  a  steady-state  reaction 
zone  followed  by  an  unsteady-state  reaction  zone. 

We  decided  to  use  histories  4,  8,  7,  and  6  to  construct  an  idealized 
set  of  particle  velocity  histories  for  the  Lagrange  analysis.  We  first 
averaged  the  noise  at  the  wave  front  on  histories  4,  7,  and  8  to  make 
the  particle  velocity  smooth  and  thereby  standardize  the  shape  of  the 
reaction  zone  at  the  wave  front.  This  smooth  peak  was  then  added  to 
history  6  to  generate  a  modified  history  6  for  the  Lagrange  analysis. 
Histories  8,  7,  and  the  modified  history  6  were  then  smoothed  to  obtain 
the  idealized  set  of  particle  velocity  histories  for  detonating  coarse¬ 
grained  Amatex  20  shown  in  Figure  9.  The  Lagrange  analysis  based  on 
these  particle  velocity  histories  was  performed  with  the  GUINSY  code^ 
written  by  L.  Seaman  of  SRI  International. 

Results 

The  results  of  this  Lagrange  analysis,  i.e.,  the  corresponding 
pressure-histories  and  the  Lagrange  compression  and  release  paths  in  the 
(p-u)  and  (p-v)  planes  at  the  three  gage  locations,  are  shown  in  Figures 
10,  11,  and  12. 

Discussion 

The  Detonation  Process  in  Coarse-Grained  Amatex  20 

The  fact  that  particle  velocity  histories  recorded  by  different 
gages  agree  for  a  period  of  time  after  the  arrival  of  the  wave  front  was 
taken  as  evidence  that  the  detonation  wave  in  coarse-grained  Amatex  20 
is  self-sustaining.  Other  properties  of  this  self-sustaining  detonation 
wave  are  obtained  from  the  particle  velocity  histories  and  the  results 
of  the  Lagrange  analysis  based  on  these  histories. 

Figures  11  and  12  show  that  the  calculated  loading  and  unloading 
paths  at  three  Lagrange  positions  lie  on  a  single  curve.  A  particle  of 
coarse-grained  Amatex  20  entering  the  detonation  wave  is  compressed 
along  the  Rayleigh  line  OR  in  the  (p-u)  and  (p-v)  planes  and  is  then 


DATE  *  26-MAR-82  SHOT  1290-2  AMATEX,  PV8.7.6.SCS.  60.59.GAGES  I-  3 

JA-1290-15 

FIGURE  9  SMOOTHED  PARTICLE  VELOCITY  HISTORIES  IN  DETONATING 
COARSE-GRAINED  AMATEX  20,  SHOT  1290-2 

These  idealized  histones  are  based  on  the  records  from  gages  8,  7,  6,  and  4 
and  were  used  as  input  in  the  Lagrange  analysis. 
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JA-1 290-16 

FIGURE  10  CALCULATED  LAGRANGE  PRESSURE  HISTORIES  AT  THE  GAGE 
POSITIONS  IN  DETONATING  COARSE-GRAINED  AMATEX  20 

These  results  are  based  on  the  particle  velocity  histories  shown 
in  Figure  9. 
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DATE  =  26 -MAR -82  SHOT  1290-2  AMATEX,  PV8,7,6,SCS.  60,59, 

JA-1290-17 

FIGURE  11  CALCULATED  LAGRANGE  PRESSURE-PARTICLE  VELOCITY 
SHOCK  COMPRESSION  AND  ADIABATIC  RELEASE  PATHS 
AT  THE  GAGE  POSITIONS  IN  DETONATING  COARSE-GRAINED 
AMATEX  20 


These  results  are  based  on  the  particle  velocity  histories  shown 
in  Figure  9. 


FIGURE  12  CALCULATED  LAGRANGE  PRESSURE-SPECIFIC  VOLUME 
SHOCK  COMPRESSION  AND  ADIABATIC  RELEASE  PATHS 
AT  THE  GAGE  POSITIONS  IN  DETONATING  COARSE-GRAINED 
AMATEX  20 


These  results  are  based  on  the  particle  velocity  histories  shown 
in  Figure  9. 


rarified  along  the  release  adiabat  RCJR  .  RCJ  coincides  with  the 
Rayleigh  line  and  represents  the  steady-state  reaction  zone;  CJR' 
diverges  from  the  Rayleigh  line  and  represents  the  unsteady  zone.  The 
unsteady  zone  is  considered  to  be  a  Taylor  wave  because  it  is  repre¬ 
sented  at  different  gage  locations  by  a  single  curve  in  the  (p-u)  and 
(p-v)  planes. 

The  particle  velocity  histories  used  to  perform  the  Lagrange 
analysis  and  the  resulting  (p-u)  and  (p-v)  paths  were  used  to  locate  the 
CJ  point  at  the  end  of  the  reaction  zone.  The  resulting  values  of  the 
CJ  particle  velocity  and  reaction  time  determined  from  the  particle 
velocity  records  are  u^j  =  1.67  mm/ps  and  ~  0.63  ps.  The  values  of 
D  =  6.85  mm/ ps  and  =  1.594  g/cm^  were  then  used  with  the  Rankine- 
Hugoniot  jump  conditions  to  calculate  the  following  values  for  the  other 
CJ  parameters:  p^j  =  182  kbar,  Vqj  *  0.4742  cm^/g,  and  Kqj  =  3.10. 

The  Role  of  AN  in  the  Detonation  Process  in  Amatex  20 

To  develop  a  better  understanding  of  the  role  of  AN  in  the  detona¬ 
tion  process  in  Amatex  20,  we  compared  the  present  detonation  parameters 
for  Amatex  20  containing  coarse-grained  AN  with  those  obtained  earlier 
for  Amatex  20  containing  baseline  AN.  The  earlier  results  were  re¬ 
examined  before  making  this  comparison.  This  reexamination  led  us  to 
conclude  that  the  state  previously  identified  as  the  CJ  state  has  a 
higher  particle  velocity  than  the  actual  CJ  state.  Reassignment  of  the 
CJ  particle  velocity  led  to  the  following  set  of  CJ  parameters  for  base¬ 
line  Amatex  20:  D  *  6.95  mm/ps,  u^j  *  1.94  mm/ps,  Pqj  =  217  kbar,  Vqj  ■ 
0.4470  cm^/g,  K^j  =  2.58,  and  x^j  ~  0.5  ps.  The  set  of  CJ  parameters  for 
coarse-grained  Amatex  20  was  as  follows:  D  *  6.85  mm/ps,  u^j  =  1.67 
mm/ps,  pgj  =  182  kbar,  v^j  ■  0.4742  cm^/g,  Kqj  *  3.10,  and  x^j  ~  0.63  ps. 
Comparison  of  both  sets  clearly  demonstrates  that  the  particle  size  of 
AN  influences  the  detonation  process  in  Amatex  20. 

The  method  of  calculating  detonation  parameters  developed  by 
Kamlet^-^  was  modified  and  used  to  interpret  the  difference  between 
these  sets  of  CJ  parameters.  Kamlet's  method  is  based  on  the  results  of 


BKW  detonation  calculations  and  consequently  assumes  that  an  explosive 
is  completely  decomposed  in  the  detonation  process.  Extension  of  this 
method  to  nonideal  detonation^  must  therefore  be  based  on  the  assumption 
that  an  explosive  exhibits  nonideal  behavior  because  its  detonation 
products  do  not  attain  chemical  equilibrium.  We  thus  assumed  in  treat¬ 
ing  Amatex  20  that  the  AN  decomposes  completely  but  that  the  reaction 
between  the  oxygen  from  the  AN  and  the  carbon  from  the  RDX  and  TNT  does 
not  necessarily  proceed  to  completion.  This  assumption  is  expressed  in 
terms  of  the  reaction  coordinates  introduced  earlier  in  this  report  by 

=  1  and  0  <  X  <1. 

4 

It  is  convenient  to  perform  two  limiting  calculations,  one  with 
X^  =*  1  corresponding  to  ideal  detonation,  and  the  other  with  X^  =  0 
corresponding  to  nonideal  detonation  with  complete  decomposition  of  the 
explosive  but  no  subsequent  reaction  among  the  detonation  products. 
Values  of  experimentally  determined  detonation  parameters  smaller  than 
those  calculated  with  X^  =  0  can  then  be  taken  as  strong  evidence  that 
the  decomposition  of  AN  is  incomplete  in  the  detonation  process.  This 
approach  to  nonideal  detonation  is  considered  to  be  more  realistic  than 
an  approach  based  on  TIGER  or  BKW  calculations^  because  it  is  not 
restricted  by  the  assumption  of  thermal  equilibrium  that  such  thermo¬ 
hydrodynamic  codes  use  to  make  detonation  calculations.  Use  of  the 
thermal  equilibrium  condition  and  the  chemical  equilibrium  conditions  to 
calculate  the  composition  of  the  detonation  products  in  such  codes 
automatically  restricts  their  treatment  of  nonideality  to  incomplete 
decomposition  of  the  explosive. ^ 

Kamlet's  semiempirical  method  for  calculating  detonation  parameters 

O 

of  CHNO  explosives  with  densities  greater  than  1.0  g/cm  uses  the 
equations 

pcj  -  PQ2  K<J)  ,  -  NM1/2  Q1/2  ,  K  =  15.58  (84) 

D  -  A<j>1/2(1  +  Bpo)  ,  A  *  1.01  ,  B  =  1.30  (85) 


the  conditions  X^  =  =  X3 
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where  N  is  the  number  of  moles  of  gaseous  detonation  products  per  gram 
of  explosive,  M  is  the  average  molecular  weight  of  these  gases,  Q  is  the 
standard  heat  of  reaction  per  gram,  and  pQ  is  the  initial  density. 

Values  of  N,  M,  and  Q  are  estimated  from  an  abitrary  assumption  about 
the  detonation  product's  composition,  so  that  calculations  require  no 
other  input  than  the  explosive's  elemental  composition,  loading  density, 
and  heat  of  formation.  The  H2O-CO2  arbitrary  assumption  about  the 
detonation  products  composition  and  the  heats  of  formation  for  RDX,  TNT, 
AN,  and  Amatex  20  shown  in  Table  5  were  used  to  calculate  the  values  of 
N,  M,  and  <f>  shown  in  Table  5.  The  values  of  a,  b,  c,  and  d  for  100  gm 
of  Amatex  20  considered  as  a  composite  explosive  are  a  =  1.5037,  b  3 
3.4216,  c  3  2.0692,  and  d  3  3.0978. 

We  use  the  superscript  e  to  denote  that  the  detonation  products 
attain  chemical  equilibrium,  and  we  use  the  superscript  f  to  denote  that 
the  composition  of  the  detonation  products  is  frozen  with  the  products 
from  each  component  in  chemical  equilibrium.  To  make  the  ideal  detona¬ 
tion  calculations  for  Amatex  20  (\^  =  1),  we  set  <J>  3  d>^Q  =  6-070  and  use 
Eqs.  (84)  and  (85)  to  calculate  the  CJ  pressures  and  detonation  veloc¬ 
ities  for  p  3  1.613  g/cw?  and  p  3  1.594  g/cm8.  To  make  the  limiting 
o  o 

nonideal  detonation  calculation  for  Amatex  20  (\^  =  0),  we  assume  that 
each  explosive  component  decomposes  into  its  own  equilibrium  detonation 

products,  set  4>  ■  “  0,2  Qrdx  +  0,4  ^TNT  +  0,4  QAN  “  4*865>  and  use 

Eqs.  (84)  and  (85)  to  calculate  the  CJ  pressures  and  detonation  veloc¬ 
ities  for  p  3  1.613  g/cm?  and  p  3  1.594  g/cm^.  The  results  of  these 
o  o~ 

calculations  are  shown  in  Table  6. 

Experimental  end  calculated  values  of  CJ  parameters  for  Amatex  20 
can  now  be  compared  to  investigate  the  role  of  AN  in  the  detonation 
process.  The  relationship  between  the  experiment  CJ  parameters  for 
baseline  Amatex  20,  (D  3  6.95  mm/ps,  pCJ  3  217  kbar,  u^j  3  1.9  mm/ps}, 
and  those  calculated  for  the  same  initial  density,  D  3  6.90  mm/ps, 

P£j  3  197  kbar,  u^j  3  1.77  mm/  ps,  3  0,  is  consistent  with  the 
following  conclusions: 


Table  5 


VALUES  OF  DETONATION  PARAMETERS  FOR  RDX,  TNT.  AN,  AND 

AMATEX  20 


Explosive 


(kcal/loOg) 


(moles/g) 


M 

/mole 


RDX,  C3H6N606 

6.62 

0.0338 

27. 

TNT,  C7H5N306 

-7.84 

0.0253 

28. 

AN,  H4N203 

-109.12 

0.0437 

22. 

Amatex  20,  C^H^N 

T3 

O 

O 

-45.5 

0.344 

26. 

Table  6 

CALCULATED  CJ  PARAMETERS  FOR 

AMATEX 

20 

V 

PCJ 

D 

Amatex  20 

(g/cnr) 

(kbar) 

(nun/ us ) 

i 

(1)  In  the  planar  detonation  process  in  baseline  Amatex  20, 

AN  with  a  mean  particle  size  of  ~0.7  mm  decomposes 
completely  in  the  reaction  zone,  but  the  detonation 
products  produced  by  RDX,  TNT,  and  AN  do  not  attain 
chemical  equilibrium. 

(2)  In  the  reaction  zone  of  the  planar  detonation  wave  in 
baseline  Amatex  20,  the  rate  of  decomposition  of  AN  is 
faster  than  the  rate  of  recombination  of  the  oxygen 
produced  by  the  AN  and  the  carbon  produced  by  RDX  and 
TNT . 

Comparison  of  the  experimental  CJ  parameters  for  coarse-grained  Amatex 
20,  {d  =  6.85  mm/ps,  p^j  =  182  kbar,  u^  =  1.67  mm/ps},  with  those 
calculated  for  the  same  initial  density,  {D  =  7.64  mm/ps,  p  =  240 
kbar,  ucj  =  1.97  mm/  ps ,  =  1 } ,  and  {d  =  6.85  mm/ps,  pCj  =  192  kbar, 

u^j  =  1.76  mm/  ps,  =  0}  does  not  lead  directly  to  a  conclusion 
because  the  experimental  values  are  approximately  equal  to  those 
calculated  for  no  reaction  among  the  detonation  products. 

A  conclusion  can  be  reached,  however,  by  assuming  that  the  reaction 
processes  in  coarse-grained  Amatex  20  are  similar  to  those  in  baseline 
Amatex  20  and  taking  account  of  the  following  facts: 

•  In  baseline  Amatex  20,  the  recombination  reaction 
between  oxygen  and  carbon  proceeds  in  the  reaction  zone 
as  the  AN  decomposes. 

•  The  reaction  time  in  coarse-grained  Amatex  20  is  longer 
than  the  reaction  time  in  baseline  Amatex  20. 

•  The  experimentally  determined  value  of  Kqj  =  3.10  for 
coarse-grained  Amatex  20  is  larger  than  the  values  of 
KCJ  resulting  from  the  detonation  calculations. 

It  follows  from  our  assumption  about  the  reaction  processes  in  Amatex  20 
that  the  recombination  reaction  between  oxygen  and  carbon  proceeds  in 
the  reaction  zone  of  coarse-grained  Amatex  20  for  a  longer  time  than  it 
proceeds  in  the  reaction  zone  of  baseline  Amatex  20.  Because  this 
recombination  reaction  liberates  more  energy  than  the  decomposition  of 
AN,  we  propose  the  following  conclusion,  compatible  with  our  first  and 
second  conclusions: 


(3)  In  the  planar  detonation  process  in  coarse-grained 

Amatex  20,  AN  with  a  mean  particle  use  of  ~2  mm  does  not 
decompose  completely  in  the  reaction  zone,  and  the 
detonation  products  produced  by  RDX,  TNT,  and  AN  do  not 
attain  chemical  equilibrium. 

The  higher  value  of  K^j  also  supports  the  conclusion  that  the  decom¬ 
position  of  AN  is  incomplete  in  the  reaction  zone  in  coarse-grained 
Amatex  20. 

Calculations  with  the  burning  rate  of  AN  written  as  a  function  of 
pressure  using  the  reaction  times  for  baseline  and  coarse-grained  Amatex 
20  are  recommended  to  estimate  the  amount  of  unreacted  AN  in  the 
reaction  zone  of  coarse-grained  Amatex  20. 
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